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■ Abstract 

■ A general framework for a systematic quasi-localization of canonical general 
^ I relativity and a new ingredient, the requirement of the gauge invariance of the 

boundary terms appearing in the calculation of Poisson brackets, are given. As a 
consequence of this it is shown, in particular, that the generator vector fields (built 

', from the lapse and shift) of the quasi-local quantities must be divergence free with 

ZJ ' respect to a Sen-type connection; and the volume form induced from the spatial 

metric on the boundary surface must be fixed. 

1 Introduction 

> ■ 

Q^ I Conserved quantities have always had distinguished role in physics. Though for a general 

^ ■ system no systematic way of finding them is known, but for systems whose dynamics 

(-r^ ■ can be described by a Hamiltonian in the canonical framework there is a way. The first 

^ ! systematic investigation of Einstein's general relativity (GR) in its canonical form was 

O I done in the ADM variables and was focused on asymptotically flat conflgurations [1]. 

^ I One of the key objects in the canonical formulation of the vacuum general relativity is 

■ the constraint function ('parameterized' by a function N and a vector field N°- on the 

■ manifold S, called the lapse and the shift, respectively): 

^['^•^1 2^(«-2^ + w[;r^''^-''-''°'])^^+ 

+ {2D,p^'')h,N'^y-x. (1.1) 

Here the canonical variables are the fields hab and p"'' on the connected n-manifold S, hab 
being the (negative definite) spatial metric, is the corresponding Levi-Civita covariant 
derivative, R is its curvature scalar and k := SnG with Newton's gravitational constant 
G, and we allow a nontrivial cosmological constant A to be present. (Though primarily 
we are interested in the physical (3+1) dimensional case, the analysis can be done in 
(n -|- 1) dimensions without any extra effort, but n > 2.) 



*Talk given at the Conference on Recent Results in Mathematical Relativity, The Erwin Schrodinger 
Institute, Vienna, 2008 August 20-21, and dedicated to Bobby Beig on the occasion of his 60th birthday. 
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One of the basic observations of Arnowitt, Deser and Misner is that the constraint 
functions play a double role in the dynamics of GR (see also [2]). In fact, the constraint 
part of Einstein's equations is equivalent to C[N^ N^] = for every and N'^ and their 
formal variational derivatives with respect to the canonical variables (given explicitly by 
fl2.2p - (l2.3l) ) appear in the canonical equations of motion, which are just the evolution 
parts of the field equations: 

where now H[N,N^] = C[N,N^]. Thus, apparently, it is the constraints that generate 
the evolution of the states of the theory through the canonical equations of motion in the 
phase space, i.e. the constraint functions appear to play the role of the Hamiltonian in the 
canonical formulation of general relativity. Another important observation of Arnowitt, 
Deser and Misner is that by the integral of the Landau-Lifshitz pseudotensor in an 
asymptotically Cartesian coordinate system, the total energy and linear momentum can 
be introduced, and these quantities turned out to be conserved during the time evolution 
of the system. 

However, as Regge and Teitelboim pointed out |3], the constraint functions C[A^, A^^] 
with the 1/r and fall-off for the canonical variables are not functionally differentiable 
in the strict sense (see e.g. [1]). The total variation of C[N, N^] with respect to hab and p""^ 
yields not only the expected volume terms, i.e. the formal variational derivatives (12. 2^ - 
(12. 3p contracted with 6hab and 6p"''', respectively, but integrals on the boundary of S at 
infinity as well (see equation (12.11) ). Thus, strictly speaking, the functional derivatives of 
C[N, N^] are the sums of smooth fields and distributions concentrated on the boundary 
of E. Therefore, if we want to recover the correct evolution equation for the smooth 
tensor fields as the Hamiltonian equations of motion (rather than some distributional 
generalization of them), then the Hamiltonian must be functionally differentiable with 
respect to the canonical variables. Since adding a boundary integral to C[N, N'^] does not 
change the formal functional derivatives, Regge and Teitelboim searched for the correct 
Hamiltonian in the form 

H[N,N^] = C[N,N']+- (f B{N,N^)dS, (1.3) 

where the integral on S is understood as the r — > oo limit of the integrals on large 
spheres of radius r in the asymptotically fiat ends of E and B{N, N'^) is some expression 
of the canonical variables and a linear expression of A^ and A^*^. They showed that there 
is, indeed, a boundary term which makes (II. 3p differentiable. Moreover, as a bonus, 
for appropriately chosen (A^, A^'^) this Hamiltonian automatically reproduces the total 
energy and linear momentum of Arnowitt, Deser and Misner as its value on the constraint 
surface and, in addition, the spatial angular momentum and centre-of-mass could also 
be introduced. Thus, the claim of a pure mathematical consistency yielded a physically 
highly desirable result. 

Nevertheless, as Beig and O Murchadha [5] showed, the asymptotic form of A^ and 
A^"^ should not be prescribed by hand. That is a consequence of the requirement of the 
compatibility of the boundary conditions and the evolution equations: Since the evolution 
equations involve the lapse and the shift, and they must preserve the boundary conditions 
imposed on the canonical variables, we get a restriction on the asymptotic form of A^ and 
A^^. In fact, in the leading order they depend on the asymptotically Cartesian spatial 
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coordinates just like the Killing vectors of Minkowski spacetime. (On the other hand, 
the time dependence of and A^^ remained unrestricted.) Without this compatibility 
condition the initial data set, satisfying the boundary conditions at the initial instant, 
would evolve in the next instant into a new data set that would not belong to the actual 
phase space, i.e. the evolution would take a data set out of the phase space. Beig and 
O Murchadha refined the Hamiltonian of Regge and Teitelboim, too, such that the new 
Hamiltonian is not only functionally differentiable, but finite on the whole phase space 
(rather than only on the constraint surface), and that these Hamiltonians form a Lie 
algebra with respect to the Poisson bracket as the Lie product. The constraints form a 
Lie ideal in this algebra, and their quotient, the algebra of observables, is isomorphic to 
the Poincare algebra. The conserved quantities then become coordinates in this quotient 
algebra. 

Since the evolution equations specify only how the lapse and the shift depend asymp- 
totically on the spatial coordinates but not on the time coordinate, moreover the spatial 
angular momentum and the centre-of-mass of Beig and O Murchadha do not transform 
in the correct way under an asymptotic Poincare transformation in the spacetime, further 
refinements of the previous results were needed. As a resolution of these difficulties in 
[HI [7] a distinction between the evolution vector fields that should be used in the Hamilto- 
nian to generate the evolution of the states and the vector fields built from the lapse and 
the shift that should be used to define ADM type conserved quantities was made. The 
latter must be the asymptotic spacetime Killing fields, which turned out to be a special 
case of the former. 

In the present paper, these ideas will be applied to the case in which the manifold 
E is compact with a smooth boundary S := 9S, i.e. at the quasi-local level (see e.g. 
\8\). Thus, now we are interested in the Hamiltonian dynamics of general subsystems 
of the universe. The extension of the investigations of canonical general relativity to 
the quasi-local case is required by solving several problems. First, to have a deeper 
understanding of the (geometrical or thermodynamical) properties of black holes, for 
example to formulate the various (geometric) inequalities or the entropy bounds for black 
holes, the conserved quantities or, more generally, the observables of the gravitational 
'field' must be introduced quasi-locally. A further motivation of searching for quasi-local 
observables is the remarkable result of Torre [9| that all the global observables for the 
vacuum gravitational field in a closed universe, built as spatial integrals of local functions 
of the initial data and their finitely many derivatives, are necessarily vanishing. Thus 
in closed universes we can associate non-trivial, locally constructible observables only 
with subsystems, bounded at one instant by some closed spacelike 2-surface. Another 
motivation is the claim to see the content of the basic existence and uniqueness results of 
Friedrich and Nagy [lO] for the initial-boundary value problem for the vacuum Einstein 
equations from the Hamiltonian point of view. 

Though the first few steps towards the systematic quasi-localization of canonical gen- 
eral relativity using the ADM variables have already been done [TT], the project remained 
incomplete and essential new ideas were needed. In the present paper, we continue these 
investigations. We refine our previous framework, and, in addition to the ideas intro- 
duced in the asymptotically flat context (and discussed above) , only flrst principles (such 
as gauge invariance in every sense, covariance, etc), but no ad hoc ideas or elements (e.g. 
some a priori reference conflguration, gauge choice or even implicitly given background 
structure) will be used. We think that these would contradict the principle of equivalence 
(see also [l2]) and hence the very spirit of Einstein's general relativity. It is the theory 
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itself that should tell us what the boundary conditions, the observables, etc are, and we 
must read off these from the structure of the theory itself. For a different view (but an 
absolutely legitimate strategy), see e.g. [T3] . 

Though in the present paper we cannot complete the quasi-localization programme, 
we raise a new issue that should be discussed in connection with the quasi-localization of 
canonical GR. Namely, we argue that the boundary terms appearing in the calculation 
of the Poisson bracket of two Hamiltonians must be gauge invariant in every sense. We 
show that the requirement of this gauge invariance yields a restriction on the lapse and 
shift: the generator vector field built from them according to K"^ = Nt"" + N"-, where 

is the future pointing unit timelike normal to the spacelike hypersurfaces S in the 
spacetime, must be divergence free with respect to the Sen-type connection induced on 
the boundary, AaK'^ = 0. Remarkably enough, this is one of the ten (or, in {n + 1) 
dimensions, the |(n + l)(n + 2)) spacetime Killing equations. A consequence of this 
condition is that the induced volume form eei...e„_i on the boundary S remains constant 
during the evolution, and hence it is natural to impose feei...e„_i = 0. feei...e„_i = 0, as 
a part of the boundary conditions for the canonical variables, has already been found in 
special contexts, e.g. in connection with the functional differentiability of the constraint 
functions. The present investigations show that feei...e„_i = should be a part of the 
'ultimate' boundary conditions, too. 

In the following section we formulate (and refine the previous attempts of) the quasi- 
localization programme of canonical GR and discuss the tools and the technical details 
that we need as well as the partial results. In section 3 we quasi-localize the canonical for- 
mulation of a single, real scalar field, and we will see that the boundary terms appearing in 
the Poisson bracket of two Hamiltonians should be interpreted as the energy-momentum 
and angular momentum flux, and hence a gauge invariant observable. In subsection 14.11 
we return to general relativity and study the analogous boundary terms in the Poisson 
bracket of two constraints. We derive here the boundary conditions AaK"" = and 
feei...e„_i = 0. Finally, in subsection 14.21 we find arguments both in favour of and against 
a Hamiltonian boundary term which is the basis of several suggestions for the quasi-local 
energy-momentum. Though there are a number of key issues even in connection with the 
present partial results (e.g. how they are related to the maximal dissipative boundary 
condition of Friedrich and Nagy [ID]), we should leave these to a future investigation. 

Our notations and conventions are essentially those that were used in [H]. In par- 
ticular, we use the abstract index formalism, the signature of the spacetime metric is 

1 — n, and the curvature is defined by —R'^hcd^'^ '■= (-Dc-Dd — DdDc)X'^. The analysis is 
based on certain formulae given explicitly in [6l[lT]. Our standard reference in canonical 
formalism is [I4j . 

2 The quasi-localization programme 

2.1 The general programme 
2.1.1 The basic requirements 

In the quasi-local canonical formulation of general relativity, the basic object is the quasi- 
local configuration space Q(S) over the connected n-manifold S with smooth boundary 
S := (9S. This is the space of smooth (negative definite) metrics hab on S satisfying certain 
not-yet-specified boundary conditions on S. The quasi-local phase space is defined to be 
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its 'cotangent bundle' T*Q(E), endowed with the natural symplectic structure. Thus 
the elements of the quasi-local phase space are the pairs (/ia6,p"*) of fields, where 
is a symmetric contravariant tensor density of weight one. (p"* is usually interpreted 
as a 1-form at the point hab € Q(S): if habiu) is any smooth 1-parameter family of 
metrics such that /lab(O) = /iafe, i-e. hab{u) is a 'curve' in through hab, then the 

tensor field 5hab on S, defined pointwise as 5hab{p) '■= iAhabiPi'^) l^'^)\u=Ot Vp G S, 
is called the tangent of the 'curve' at hab- In fact, the directional derivative of any 
functionally differentiable function F : M along the 'curve' at the point hab is 

5F := {(\.F{hab{u))/(\u)\u=Q = Jj.-^f-^dhabd'^x, which defines the action of 5hab on any 
such F and the natural pairing of 6hab and the 'exact 1-form' dF = 4^ at the same 
time. Thus, it is natural to define the action of p""^ on the 'tangent vector' 6hab at the 
point hab £ by the integral 5habp"''^d^x. Note that while in the asymptotically 

flat context the requirement of the finiteness of the integral J^, Shabp°'^d'^x together with 
the evolution equations restrict the fall-off rate k of the metric, namely [6] it must be 
k > |(^^ — 1), in the quasi-local case we obtain no restriction for the canonical variables.) 
Clearly, the differentiability of functions on the quasi-local phase space depends not only 
on the function itself, but on the boundary conditions that are imposed on the canonical 
variables in the definition of T*Q(E). We stress that the boundary conditions are parts 
of the definition of the phase space. (For a more detailed discussion of this issue, see e.g. 

The canonical symplectic structure can be characterized equivalently by specifying 
the Poisson bracket of any two functionally differentiable functions. For any such G and 
H : T*Q(S) -> M, it is {G, H} := - f.£;)d"x. 

In a more pedagogical approach the lapse and the shift are also considered to be 
configuration variables (see e.g. [IE]), and there are additional constraints (the 'primary 
constraints') that the momenta conjugate to the lapse and the shift are vanishing. Then 
a systematic constraint analysis shows that this 'big' phase space T*Q(S) can always be 
reduced to T*Q(S) in a straightforward way, and the role of the lapse and the shift is 
reduced from dynamical variables only to 'parameters'. In this sense, T*Q(S) can also 
be considered as a 'partially reduced' phase space. We save this reduction by starting 
with T*Q(S) as the phase space. 

Our ultimate aim is the quasi-localization of canonical general relativity, i.e. finding 

1. a boundary term B{N, N^) (we will call its integral the Hamiltonian boundary term), 

built from the canonical variables and depending linearly on the lapse and the shift; 

2. boundary conditions for the canonical variables {hab^P""^) on iS; 

3. boundary conditions for the lapse and the shift on S 

such that 

i. the Hamiltonian H : T*Q(S) — M, given by (11.31) and 'parameterized' by lapses and 

shifts satisfying the boundary conditions in point 3, is functionally differentiable 
with respect to the canonical variables; 

ii. the evolution equations (11.21) with lapses and shifts satisfying the boundary condi- 

tion in point 3. above preserve the boundary conditions imposed on the canonical 
variables in point 2; 
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iii. the constraints close to a Poisson algebra C (which we call the quasi-local constraint 

algebra); 

iv. the value of the Hamiltonian on the constraint surface (i.e. if C[N, N^] = 0) must 

be a 2 + (ra — l)-covariant, gauge-invariant expression of the boundary data on S. 

Before turning to the mathematical realization of these requirements we should discuss 
these points. 

2.1.2 The discussion of the requirements 

Clearly, requirement (i) is just that of Regge and Teitelboim [5, and (ii) is taken from Beig 
and O Murchadha While the first is a generally accepted condition in the relativity 
community, the second is apparently not appreciated enough. Moreover, the usual (and 
in the quasi-local canonical approaches almost exclusively adopted) view is that the lapse 
and the shift on the boundary S should be arbitrary, just because they are thought of 
as the n + 1 pieces of the spacetime evolution vector field, and their arbitrariness are 
thought to express our freedom to choose the fleet of observers in the spacetime as we 
wish. 

However, as we learnt from the structure of canonical GR of asymptotically flat space- 
times, requirement (ii) may yield non-trivial restrictions on the lapse and the shift. More- 
over, as the study of the quasi-local conserved quantities of matter fields in Minkowski 
spacetime shows (and as we summarize the basic things in subsection 13. II) . we must make 
a distinction between the symmetries (which define the conserved quantities) and the 
general evolution vector fields (which are used to define the time evolution of the initial 
data set in the spacetime). The former is the solutions of a linear partial differential 
equation (and the number of their independent solutions is finite), but the latter is quite 
arbitrary. Thus, in general relativity, the lapse and shift parts of a general evolution 
vector field may be, but the lapse and shift parts of a vector field defining 'conserved' 
quantities are probably restricted. In particular, as we will see in subsection 12.3.31 the 
'most natural' quasi-local Hamiltonian cannot provide an acceptable expression for the 
energy-momentum if the lapse and shift on S are chosen independently of the canonical 
variables. We expect that these conditions are given only implicitly, possibly in the form 
of a system of linear partial differential equations (p.d.e.) on S. In fact, in section H] 
we derive an equation for the lapse and shift which should be a part of such a system of 
linear p.d.e. 

Requirement (iii) expresses the claim that the gauge content of the theory must be 
represented in a correct way by the constraints even quasi-locally. The motivation of this 
expectation comes from the applicability of the standard canonical formalism (see e.g. 
[T4]). In fact, for n > 3 general relativity is expected to have internal degrees of freedom, 
which are expected to be represented in the canonical framework by the points of the 
reduced phase space. However, strictly speaking the reduced phase space can be formed 
as the space of gauge orbits in the constraint surface only when the Hamiltonian vector 
fields of the constraint functions form an involutive distribution or, equivalently, if the 
constraint functions close to a Poisson algebra. 

Though in the asymptotically flat case the Hamiltonians of Beig and O Murchadha 
also close to a Poisson algebra (in which the constraints form an ideal) [5], we will see 
in subsection 13.2.21 that quasi-locally this is not true even for the scalar field, and hence 
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it cannot be expected in general relativity either. In fact, the boundary terms appear- 
ing in the Poisson bracket of two Hamiltonians (and which destroy the Poisson algebra 
structure) describe incoming and outgoing energy-momentum and angular momentum 
fluxes. Thus quasi-locally the algebra of observables is much bigger than the set of the 
Hamiltonians. 

In the present investigations, we do not use any 'natural' (or rather ad hoc) conditions 
e.g. for the canonical variables. We use only first principles, such as gauge invariance (in 
any sense) or covariance. Requirement (iv) is the manifestation of this idea. In particular, 
the boundary term B{N, N'^) should not depend on the normal directional derivatives of 
the lapse and the shift, or on the actual lapse and shift separately, but only on their 
spacetime-covariant combination K'^ := Nf^ + N'^. Moreover, it should not depend on 
the actual choice for the normals (f^jf") of S but only on the intrinsic and extrinsic 
geometry oi S (but not e.g. on the extension of the geometrical quantities off 5). 

However, because of the strong interplay between the Hamiltonian boundary term 
and the boundary conditions both for the canonical variables and the lapse and shift, it 
would be extremely difficult to determine these three unknown things in a single step. 
Indeed, in the traditional approach usually one of these three is a priori given, and by the 
requirements (i) and (ii) we can determine the other two in a relatively straightforward 
way. For example, in the asymptotically fiat case we have a priori boundary conditions, 
from which Beig and O Murchadha could determine both the Hamiltonian boundary term 
and the boundary conditions for the lapse and shift. In subsection 12.31 we summarize the 
special cases in which various Hamiltonian boundary terms were a priori given, from 
which the boundary conditions could be determined. (The necessary tools and technical 
ingredients will be summarized in subsection 12. 2[ ) Thus, in order to be able to manage 
the general problem, a new idea would be needed. 

As we already mentioned in the introduction, such a new idea could be the observation 
that the boundary terms in the Poisson bracket of two Hamiltonians (which we call the 
Poisson boundary term) should be the sum of the Hamiltonian boundary term and terms 
analogous to the energy-momentum and angular momentum fiows between the system 
and the rest of the universe, and hence these must be gauge invariant. This idea is based 
on the moral of the analogous investigations of a scalar field, the subject of section 3. 

2.2 The main ingredients 
2.2.1 The variational formula 

Though formally the manifold S on which the canonical variables {hab,P"'^) as fields are 
defined is an abstract n-manifold, it is convenient (and illuminating, and hence useful) 
to consider this as the typical leaf of a foliation {Sf} of (a piece of) the spacetime by 
spacelike hypersurfaces and identify S at the coordinate time instant t with E^. Thus 
although in the phase-space context there is no reason to speak about the normal of 
S, in the spacetime we have a uniquely defined future pointing unit timelike normal 
of the leaves (and hence a projection P^" := 6^ — f^tb of the spacetime tangent 
spaces to the tangent spaces of the leaves, too), and we frequently rewrite formulae given 
in the phase-space context to its spacetime form or vice versa freely if needed. For 
example, the spacetime form of the lapse and the shift is the so-called evolution vector 
field = Nt°- + N"-, by means of which the spacetime form of the constraint function 
(II. ip is just the integral C[^^] := \ jj^i°'{^^Gab + A(y'afe)t^dS. Here Qab is the spacetime 
metric, ^^Gab is the corresponding Einstein tensor and the induced volume element on E 
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is dS := ^t^eeai...a„- {£ai...a„+i IS the volume {n + l)-form in the spacetime.) Note that 
the timelike unit normal to S is globally well defined if the spacetime is time orientable. 

Also, we will need the expression of the canonical momentum in terms of the Lagrange 
variables, i.e. the configuration and velocity variables. The latter is defined with respect 
to some spacetime vector field such that the hypersurfaces of the foliation are 
obtained from S by Lie dragging along the integral curves of We assume that these 
hypersurfaces are spacelike. Let t"" be their future pointing unit timelike normal, the 
lapse and the shift parts of ^"^ are := te^^ and A^*^ := PfC,^, respectively, and the 
acceleration of the hypersurfaces is Oe := fVate = —DglnN. Here denotes the 
induced Levi-Civita derivative operator on TS^. The time derivative of a purely spatial 
tensor field T^[" is defined to be the projection to the hypersurfaces of its Lie derivative 
along i.e. fl^- := {L^T;^-)P^...P^... = N{LtTl-)P^ ...P^ .. + L^T^^- . In particular, 
hab = 2Nxab + Lf^hab- Thus essentially the extrinsic curvature Xab of the leaves of the 
foliation plays the role of the velocity, by means of which the canonical momentum is 
known to be p"'^ = ^^/ifAix""^ ~ X^"^)- Here x is the /iab-trace of Xab- 

To calculate the total variation of (jl.ip . let A^(m), N"-{u), hab{u) andp"''(u), u G (— e, e), 
be any smooth 1-parameter families of lapses, shifts, metrics and canonical momenta, 
respectively. Then we define the corresponding variation of any of their function, F = 
F{N,N'',hab,p''''), to be its w-derivative at m = 0, i.e. 6F := ldF{N{u),N''{u),hab{u), 
p°'''{u))/du)\u=o- The corresponding variation of the constraint function C[A^, A^*^], taken 
from is 



Jt, ^ ^i^ab Op^ / 

^ / I Nih'^'v^iDMab) - v^D'Shab)) + {v^D^N - h''%'D,N)6hab 



2k 



+ 2 



Here is the outward pointing unit normal of iS in S, diS is the induced volume element 
on the boundary, and the formal variational derivatives are 



{2N''VeP''^ - N^VeP''^)8hab + AKNaV^,-^ \dS. (2.1) 



ohab 2k IV \h\ ^ 

—Kap'^'^p^^)) + D^D^N - h^^D^D^N] - L^p"^ + 

n — 1 ' J J 

+ -Nh^'y^^R -2X + -jj^i^/ - rtc.)), (2.2) 

^-^^^kr^ ■■ = -%N(p,b - -^p-'K^hab) + L^Kb. (2.3) 
op ^J\h\ V n-l J 

Here Rab is the Ricci tensor of De, and note that the last line of (12. 2p is — times the 
integrand of the constraint function C[N, 0]. 

Therefore, C[N^ N^] is functionally differentiable (in the strict sense of [4]) with respect 
to the canonical variables only if the boundary integral in (12. ip is vanishing, in which case 
the functional derivatives themselves are given by (12.21) and (12. 3p . Then, provided the 
constraints are satisfied, the vacuum evolution equations (with the cosmological constant) 
are precisely the canonical equations of motion ( 11.20 . 
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2.2.2 A quick review of the geometry of the boundary surface 



To evaluate the boundary terms in (12.11) . it seems useful to split the variation of the 
metric hab at the points of S with respect to the boundary. Moreover, in the subsequent 
subsections several expressions on S, obtained originally in the traditional n + 1 form, 
will have to be rewritten in a 2 + — 1) form. Nevertheless, most of these notions have 
a non-trivial meaning only if n > 3; and hence when we use them we assume that n > 3. 
Thus now, in a nutshell, we summarize the basic geometric objects that we need in what 
follows. A more detailed discussion of these concepts is given e.g. in [T7l [8]. 

To avid confusion, the Kronecker delta on S will be denoted by P^, and the hab- 
orthogonal projection to S is := + v"'Vb. Then the induced metric and the cor- 
responding intrinsic Levi-Civita covariant derivative on S is qab '■= hcdJ^'a^b ^^'^ "^e, re- 
spectively, and let us introduce another derivative operator simply by Ag := lilDj. The 
extrinsic curvature of 5 in S will be defined by Vab '■= ^a^b^c^d- The difference of these 
two derivative operators is just the extrinsic curvature: AeA" = SeiJl^X'^) — v"'5e{ybX^) + 
(z/gfef — i'e"'Vb)X'^ for any X°' tangent to S. The induced volume — 1) form and volume 
element on S are £ei...e„_i := i''^^^£^afeei...e„_i and dS := jj^t''v''eabei...e„.i, respectively. 
Note that with these conventions, the Gauss theorem takes the form DaX'^dT, = 
— §gVaX"'dS for any vector field A" on E. 

The boundary S = dT, can be considered as a submanifold in the spacetime, too. In 
the spacetime context the induced metric is qab = gcdX^'a^bi the area 2-form on the 
2-planes normal to S is -^Eab '■= ta^b — h^a- Here, both the projection 11^ and the area 
2-form -^Sab are independent of the actual choice of the normals Note that the 

normals are not specified by 5, but if S is considered to be the boundary of S, then they 
are chosen as in the previous paragraph. If S is orientable and at least a neighbourhood 
of 5 in M is space and time orientable, then {t"-, f ") can be chosen to be globally defined, 
yielding a global trivialization of the normal bundle NS of S in M. The two derivative 
operators Ae and 5e, acting on any Lorentzian n + 1 vector field A", are defined by 
AeA'* := n^V/A'^ and ^A" := n;jAe(n^A^) + {5^ - n^)Ae((5,^ - K)X''). The extrinsic 
curvature tensor of 5 in M is Q°'eb '■= — n^AgllJ; = r^-eh — J^"'eVb, where Tab and Uab 
are the individual (symmetric) extrinsic curvatures of iS in M corresponding to the unit 
normals ta and Va, respectively. The corresponding traces are r := Tabq"''' and z/ := Uabl""^, 
respectively. The difference of the two derivative operators is this extrinsic curvature 
tensor: AgA"^ = (5eA" + {Q'^eb — Qbe"')X^ . 5^ is not only the Levi-Civita derivative 
operator 5e on the tangent bundle of S, but it acts on the normal bundle NS of 5, 
spanned by the two normals t"- and f", as well. Its action can be characterized by the 
connection 1-form Ae := {^etajv"" = {5eta)v"'. On the other hand, for vectors tangent to 
S in M the two derivative operators Ae and Ae coincide, which fact will be used several 
times when we rewrite expressions given in the n + 1 form into its 2 + (ra — 1) form. 

At the points of 5, the splitting hab = qab — VaVb implies the variation 5hab = Sqab — 
VaSvb — VbSva- It is Straightforward to determine the various projections of 6hab (for details 
see [H]). These are 



6hedK^t = Sq^dK^t, SKdv'Ut = -Sv^qab, 

6hcdV%'^ = 2v''6va = -2va5v\ (2.4) 

Thus the independent variations can be represented by dqcd^a^f 

The curvature of the connections 5e and Ae, respectively, are given by 
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Fbcd = -^e\{5,Ad - 5dA,) + ^R\,d, (2.5) 

F'^bcd = f°'bcd — Sc{Q"'db — Qbd"') + Sd{Q"'cb ~ Qbc"') + 

+Q'^ceQbd'^ + Qec'Q^db — Q"' deQbc — QedQ" do-, (2.6) 

where '^R°'bcd is the curvature tensor of the intrinsic Levi-Civita connection Se of {S,qab)- 
Of course, for n = 3, it can also be written in the form ^^Rijl'^qbd — n^g^c), where '^R is 
the curvature scalar. The curvature F°-hcd turns out to be just the pull-back to S of the 
spacetime curvature 2-form: -F'^bcd = R"'bef^l^d- various projections, 



^a^b-Fefcd — Rabcd + TacUd — TadUc — ^ac^bd + ^ad^bc, (2.7) 

t^-UlFafcd = ScTdb - SdTcb + AcVdb " AdVcb, (2.8) 

V^IllFafcd = Scl^db - Sdl^cb + A^Tdb - AdTcb, (2.9) 

fv^Fabcd = Teciy'd " Tedl^'c + ^cAd - ^d^,, (2.10) 

are the so-called Gauss, Codazzi-Mainardi, and Ricci equations. 



2.3 Special cases 

2.3.1 The quasi- local constraint algebra 

A special, genuinely quasi-local case in which the programme could be completed is when 
there is no Hamiltonian boundary term, i.e. when we are interested in the boundary 
conditions both for the canonical variables and for the lapse and shift that make the 
constraint functions functionally differentiable and close to a Poisson algebra. The sig- 
nificance of this special case is that the constraints represent those parts of the field 
equations that are expected to generate the gauge motions in phase space. Thus to un- 
derstand the gauge content of GR at the quasi-local level, we should first clarify this 
special case. 

Decomposing the boundary terms in (12.11) with respect to the boundary S according 
to subsection 12.2.21 we can read off the condition of the functional differentiability of 
C[N, iV^] [H]: the lapse and the shift are vanishing on S and the induced volume (n — 1)- 
form £ei...e„_i is fixed on S. It is straightforward to check that the boundary condition 
fcei...e„_i = is preserved by the evolution equations with shifts and lapses vanishing on 
S. Thus the quasi-local Hamiltonian phase space T*Q(S) is split into disjoint sectors 
T*Q(S, eei,...,e„_i), labeled by the value of the volume form on S. The constraint functions 
are differentiable in the directions tangential to these sectors, but not in the directions 
transversal to them. The Poisson bracket of any two constraint functions C[N, N'^] and 
C[N,N^], 'parameterized' by lapses and shifts that are vanishing on 5, is 

^^C[N,N''],C[N,N'']^ =C[N^DeN - N^DeN,ND''N - ND^N - [N,N]'']. (2.11) 

Furthermore, the new smearing fields NW^N - NW^N and ND^N - ND^N - [N, N]" 
are also vanishing on the boundary S. Therefore, the constraints close to a Poisson 
algebra C. 
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Geometrically A^l^ = 0, A^'^l^ = correspond to an evolution vector field = 
t"-N + A^" in the spacetime that is vanishing on 5, and hence the 1-parameter family 
of diffeomorphisms (pt generated by leaves S fixed pointwise. This (j)t maps S into 
a family of Cauchy surfaces for the same globally hyperbolic domain -D(S) with the 
same boundary dT^t = S. According to Bergmann [18], the gauge-invariant content of 
general relativity is the spacetime geometry, and hence any two sets of information that 
specify the same spacetime geometry must be considered to be gauge equivalent. In 
particular, two Cauchy data sets determining the same globally hyperbolic domain are 
gauge equivalent in this sense. (For different interpretations see e.g. [El [20].) Therefore, 
the evolution vector fields that are vanishing at S are precisely the generators of gauge 
motions of the actual quasi-local state in the spacetime. 

2.3.2 The algebra of the basic Hamiltonians 

If in the quasi-local Lagrangian phase space we choose the Lagrangian L := ^ /^(-R ~ 
2 A + XabX""^ ~ X^)^VWd"x, then for the basic Hamiltonian we obtain 



Thus we already have a nontrivial a priori Hamiltonian boundary term, in which both 
VatbQc^^ and Aa depend on the actual choice for the normals (t",f") of S. This Hamil- 
tonian can be made 2 + (n — l)-covariant if is restricted to be tangent to S and, to 
cure the 5*0(1, l)-gauge dependence of the connection 1-form, if 5eK^ = is required. 

In fact [11], evaluating the boundary terms in the total variation of Hq we obtain 
that Hq is functionally differentiable if the lapse is vanishing on iS, the shift is tangential 
to S, and the volume {n — l)-form is fixed on S. However, this boundary condition 
is preserved by the evolution equations precisely when SeN*^ = is satisfied. With 
these boundary conditions the basic Hamiltonians form a Poisson algebra T^o, in which 
the quasi-local constraint algebra C is an ideal. By evaluating the basic Hamiltonian 
on the constraint surface we get a function on the algebra Tio/C of observables, which 
provides a representation of the Lie algebra of the (5e-divergence-free vector fields on S. 
Though the observable OfiV*^] := —\§g N^^A^dS behaves as angular momentum in certain 
special situations (see [HI El]), but this can be non-zero even in Minkowski spacetime. 
This shows that the boundary term of Hq should probably be present in the 'ultimate' 
Hamiltonian, but still further terms are needed. 

Geometrically, N\s = 0, VaN'^\s = correspond to evolution vector fields tangential 
to S on S. The corresponding 1-parameter family of diffeomorphisms still maps -D(S) 
onto itself and preserves the boundary S, but not pointwise. Its action on S preserves 
the volume {n — 1) form. Thus Hq is certainly not the 'ultimate' Hamiltonian, it is only 
an improved version of the constraint functions. 

2.3.3 On the difFerentiabihty of the 'improved' basic Hamiltonian 

The 'bad' gauge dependence of the basic Hamiltonian (I2.12p can be improved slightly 
by hand by adding Nu = KHaV^Qc'^'^ to the integrand of the boundary integral. (Here 
'■= Vab(t^ , the trace of the extrinsic curvature of 5 in E.) The resulting expression. 




(2.12) 
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Hi[K'] ■.= C[K'] - ^j^K'^i^^SabQc' + A^dS, (2.13) 

has been derived in different forms by several authors [22l [23l [24| and used p5] to define 
quasi-local energyQ Here the first term in the boundary integral became 2 + (n — 1) 
covariant, but the second term still depends on the 5'0(1, 1) boost gauge. To cure this 
dependence, we still must require deK"^ = 0. In subsection 14.21 we will see that the 
boundary term of f l2.13p emerges naturally among the boundary terms in the calculation 
of the Poisson bracket of two constraint functions (or Hamiltonians). Moreover, in that 
context there is a natural way of curing its 50(1,1) boost-gauge dependence. On the 
other hand, as we will see, it does not seem to represent correctly the 'composition law' 
of the lapses and shifts. 

Calculating the total variation of Hi, we can determine the condition of its functional 
differentiability. However, it is enough to calculate the total variation of the 'correction 
term' Nu and to use the expression (4.2) of |TT] for the total variation 6Hq[K^] of the 
basic Hamiltonian. The total variation of Nl'^/\q\ is 

5{Niyy/J^\) = (z/5iV + ^Nv'{DJK,)q-'' + ^uNq^'dqab + 
where we used decomposition (12. 4p . This, together with SHq[K'^], yields 

-■^ y I Z/ 5iV + (Ae - TVe)5N' + VeN^AaSv" - TVeN'vJv" + 

+ {uN - TV,N' + v'D.N + XcdV%%eN')q''''yqabYs. 

Thus if the variations of the lapse and the shift were independent of the variations of 
the canonical variables, then the differentiability of Hi[K^] with respect to the metric 
hah could be ensured by keeping the whole n-metric fixed on S: 5hab\s = 0. However, 
this condition is not 2 + (n — l)-covariant, and from the compatibility of this bound- 
ary condition with the evolution equation for hab it follows that K'^ on S must satisfy 
v'^v'^V (cKd) = 0, njjf'^V(ci^d) = and WJlfV [cKd) = 0, which are conditions on the 
derivative of the lapse and shift in the direction v"" normal to S, too. Another possibility 
is that the induced metric qab on S is fixed and the shift is tangent to S. Apparently, 
this is a weaker condition for (A^, A^") than what we had in subsection I2.3.2[ However, 
VeN^\s = is invariant with respect to the S0{1, 1) transformation of the normals (t", f") 
only if we require A^l^ = 0, too, i.e. we arrived back to the basic Hamiltonian. Hence, we 

^The same terms also appear in the general expression of the quasi-local quantities based on M0ller's 
boost-gauge invariant, but still 0(l,n) gauge-dependent superpotential, and hence in the spinorial ex- 
pressions based on the Nester-Witten 2-form in 3+1 dimensions, too. For details see the forthcoming 
updated version of [8]. 
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conclude that the Hamiltonian 112. 13\} could be the 'ultimate ' quasi-local Hamiltonian only 
if the lapse and shift are not independent of the canonical variables. Therefore, according 
to our expectation in subsection 12. 11 the lapse and the shift should satisfy a certain linear 
differential equation on S. In subsection 14. II we derive such an equation, but to motivate 
those investigations first we study the quasi-local canonical formulation of a single real 
scalar field. 

3 Illustration: Matter fields 

3.1 Conserved quantities and flux integrals for general matter 
fields 

Let the matter fields be described by the symmetric energy-momentum tensor T^b, which 
is divergence free if the field equations are satisfied. Let K'^ be any vector field, S a 
smooth compact spacelike hypersurface with smooth boundary S := dTi, and let us form 
the integral 

Q^K^] := f K^T^^^Sbe,...^. (3.1) 

Let be another, arbitrary smooth vector field on M, define to be the 1-parameter 
family of hypersurfaces by Lie dragging S along the integral curves of such that Eq = S 
and form the 1-parameter family of integrals (13. ip on these hypersurfaces. Then the 
derivative of these integrals with respect to the natural parameter t along the integral 
curves of i^*^ at t = is 

However, the integrand of the second integral on the right can be rewritten into the 
exact n-form ~ (n-i)! ^[ei (^e2...e„_i]abC^'"^-^c), and hence by the Stokes theorem it can be 
converted into a boundary integral. Thus finally we have 

Therefore, if the energy-momentum tensor is divergence free, there are two roots of the 
non-conservation of the quantity Qs [K^] : The non-Killing nature of the vector field 
and the boundary integral. 

If K"^ is a Killing field, then the vanishing of the right hand side of (13.21) can be 
expected only for tangent to S. The hypersurfaces that such a generates are 
such that the boundaries of all these coincide, dT,t = S, and they are simply other Cauchy 
surfaces for the same globally hyperbolic domain -D(S). Therefore, the Qs[-f^^] for Killing 
K'^ must in fact be conserved for tangential to S on S. For example, if (M, gab) is the 
Minkowski spacetime with Cartesian coordinates {x-}, a = 0, then the general 

Killing field has the form = VeX- + Mab{x-VeX- —x-VeX-) for some constants 
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and Mab = M[ab]. Then the coefficients of these constants in Qs[i^1 =: T^P- + MabJ-- 
define the quasi-local energy-momentum and angular momentum of the matter fields. 
These are conserved during the evolution with any ^'^ being tangent to S at S, and 
transform in the correct, expected way under Poincare transformations of the Cartesian 
coordinates. Hence the quasi-local quantities can be thought of as being associated with 
S or with the whole Cauchy development -D(S) of S. 

To understand the meaning of the boundary integral in f l3.2p . suppose that is not 
tangent to S. Since the area 2-form -^Sab annihilates the part of tangential to S, 
without loss of generality we may assume that is orthogonal to S. Let B denote the 
union of the boundaries dT,t for all \t\ < e for some positive e, i.e. the hypersurface that S 
sweeps out by Lie dragging along the integral curves of (This -B is a smooth, regular 
hypersurface only if ^'^ is nowhere vanishing on S. At the zeros of ^'^ the boundaries 9Sf 
intersect each other, and at these points B collapses to {n — 1) dimensional.) Then, by 
construction, ^"'^Eat is a (not normalized) normal 1-form to the (regular parts of) B on 
S. Thus, the integrand of the boundary integral on the right hand side of fl3.2p is the flux 
density of the current T'^'^K^ through i? at 5 weighted by the 'length' of Therefore, for 
small enough At, §gi''^eabT^^KcdS/\t is the flux of the current T"^ Ki, through B between 
(9So and dT^At- The root of the non-conservation of the quasi-local quantities Qsi-f^*^] 
even for Killing K"^ is that the actual system, surrounded at a given instant by S, is not 
closed, and there can be non-trivial incoming and outgoing flows of energy-momentum 
and angular momentum. 

In particular, if is the time translational Killing held, Ke = VeX°, then Q^f-f^^] = 
P°, the quasi-local energy. If T"'' satisfies the dominant energy condition, then this is 
non-negative (and zero if and only if T"'' is vanishing on D(S)). Then by (13.21) is 
non-negative for outward pointing (i.e. for which v^^'^ < 0) spacelike or null non- 
positive for inward pointing (i.e. Ve^'^ > 0) spacelike or null and does not have a 
definite sign for timelike C,'^. In the first case only the incoming energy flow can cross B, 
yielding energy gain; in the second only the outgoing energy flow can cross B, yielding 
energy loss; while in the case of timelike B, both incoming and outgoing energy flows 
may be present. 

As a conclusion, first, a distinction between evolution vector fields generating the 
(e.g. time) evolution of the state and the generators of the quasi-local (conserved) quan- 
tities K'^ must be made. (This view was already adopted already both for matter fields 
in [HI [26] and for gravitational fields that are asymptotically fiat at spatial infinity in 
[HIT].) Moreover, the boundary integral appearing in the 'time' derivative of the 'con- 
served' quantities describes the fiux of the incoming and outgoing energy-momentum and 
angular momentum fiows. We will see in subsection 13.2.21 that a detailed and systematic 
quasi-local Hamiltonian analysis of a single real scalar field exactly reproduces the result 
(13. 2p . where the boundary integral emerges as the boundary term in the Poisson bracket 
of two Hamiltonians. Though the gravitational 'field' does not have any well-defined 
energy-momentum density but admits a Hamiltonian formulation, an analogous result 
may be expected for the gravitational 'field' as well: the Poisson boundary term must 
have a physical meaning (and hence must be gauge invariant) for appropriately defined 
'quasi-symmetry generators' K'^ on S. 
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3.2 Quasi-local Hamiltonian description of the scalar field 



3.2.1 The quasi- local phase space and the Hamiltonian 

Let $ be a real scalar field on M, whose dynamics in the spacetime is governed by the 
Lagrangian L := |5'''''(Va$)(Vb$) — where the potential V = V{^) is a given algebraic 
function of the scalar field. For the sake of concreteness, we may assume that this has the 
form V = |m^$^ + i.e. the scalar field is of rest-mass m and /i is its self-interaction 

parameter. The covariant field equation and the energy-momentum tensor, respectively, 
are 



VaV"$ + ^ = 0, (3.3) 

Tab = (V„<f ) ( V5$) - ^gab ( Ve$) {V'<l>) + gabV. (3.4) 

Tab with the explicit form of the potential V above satisfies the dominant energy condition 
precisely when /x > 0. 

The basis of the quasi-local canonical formulation of the theory of scalar field is 
the quasi-local configuration space Q(S), the space of the smooth real scalar fields on 
E satisfying certain, not-yet-specified boundary conditions. The quasi-local Hamilto- 
nian phase space is its 'cotangent bundle' T*Q(S), endowed with its natural symplectic 
structure. The canonical momenta are scalar densities 11 on S. Using the Lagrangian 
L : TQ(S) M, defined by L := CN ^/ihld^x and considered to be the function of the 
Lagrange variables ($, $), the standard canonical formalism yields for the momenta that 
n = r{Va^)^\ = ^(^ - NW^^)^\ and for the Hamiltonian H : T*Q(S) ^ M, 
introduced by H := fj^iu^ - CNy^Dd^'x, that 

H = J^{n{\^^ - Ih'^'iDa^) (Db^) + V^y^\+ iV^(De<l>)n}d"a:. (3.5) 

A straightforward calculation shows that the coefficient of is just the energy den- 
sity part fi := Tabt"'t'^ and the coefficient of A^"^ is just the momentum density part 
Je := P^Tabt^ of the symmetric energy-momentum tensor. Thus we can also write 
H = Jj, K'^Tabt'^dT,, where K"^ = Nt"" + A^"^. If more than one Hamiltonians 'parame- 
terized' by different lapse-shift pairs are considered, then to indicate which lapse-shift is 
used we write the Hamiltonian as H[N^N^] or H[K^]. 

Since H depends on the 'parameters' A^ and A^"^, the spatial metric hab and the mo- 
mentum variable H algebraically, H is functionally differentiable with respect to them, 
independently of the boundary conditions on S. (Though in the present context the 
functional differentiability with respect to A^, A^*^ and hab does not have any significance, 
in subsection 14.1.51 where we consider the Einstein-scalar system, hab will be the gravi- 
tational configuration variable, and hence the functional differentiability with respect to 
hab will be important.) The corresponding functional derivatives themselves are 

where aab '■= TcdPaPbi the spatial stress part of the symmetric energy-momentum tensor 
O, and 
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Nevertheless, the condition of the functional differentiability with respect to $ is 

^ v^'i^Nall- N{Da<^)y<^dS = 0. (3.8) 

(Here 11 is the 'de-densitized' canonical momentum: fl =: nA/|7i|.) (13.81) can be ensured 
either by fixing the configuration variable $ on S, 6^\s = 0, or by requiring the van- 
ishing of the coefficient of 5$ in (13.81) . Under any of these conditions H is functionally 
differentiable, and its functional derivative with respect to the configuration variable is 

^ = -Da{N'^fl - iV(Z)»$) v^) + N^^\. (3.9) 
Then the canonical equations of motion are 

^ = ^ = NU + N''Da<l>, (3.10) 

n = = -Da[N^\D-<^ - iV^n) - N^^\. (3.11) 

The first is equivalent to the definition of the time derivative of $, while the second to 
the field equation ( 13.31) . 

Returning to the boundary conditions, the first, 6^\s = 0, is analogous to the Dirich- 
let boundary condition in electrostatics. However, according to our requirement (ii) in 
subsection 12.1.11 the evolution equation (I3.10p must preserve this, and hence we obtain 
that at S the canonical variables must satisfy 

NU + N^Da^ = 0. (3.12) 

Thus, the canonical momentum (weighted by the lapse) must be linked to the derivative 
of $ in the direction of the shift. In particular, for vanishing shift 11 must be vanishing 
on S. The other possible boundary condition coming from (13. 8p is the requirement of the 
vanishing of the coefficient of 5$ on iS: 

v^Nall - Nv^Da^^ = 0. (3.13) 

Thus, the normal directional derivative of $ (weighted by the lapse) is linked to the 
canonical momentum on S. This is analogous to the (generalized) Neumann-type bound- 
ary condition in electrostatics. (The Dirichlet- and Neumann-type boundary conditions 
appear naturally in a covariant phase-space context too. For details, see pTl l28|.) 

3.2.2 The Poisson boundary term and the flux integral 

Let (A^, A^^) and (iV, iV'=) be two lapse-shift pairs, assume that both H[N, N""] and H[N, N"] 
are differentiable and let us calculate the Poisson bracket of two Hamiltonians parame- 
terized by them. By integration by parts, a direct calculation yields that it is 
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^^H[N,m],H[N,N^]^ = 

= H iN^DaN - N^DaN, ND^N - ND^N - [N, N] "1 + 



(3.14) 




Js 

Note that the new lapse and shift that parameterize the Hamiltonian on the right are 
exactly those that appeared in the constraint algebra of Einstein's theory (see subsection 
12.3.11) . But in addition to the Hamiltonian H the spatial integral of the spatial stress, 
contracted with the Killing operators acting on the shift vectors (and weighted by the 
lapses), also appears. These operators can be replaced by the spacetime Killing operators 
acting on appropriately defined spacetime vector fields. Indeed, let us fix a foliation of 
the spacetime with lapse M and a compatible evolution vector field := Mf^ + M^, 
where is the future pointing unit timelike normal of the leaves of this foliation, and 
define K'' := Nf" + and K'' := Rf' + N". Then the complete n+1 decomposition of 
the Killing operator V(aA'fe) with respect to this foliation and evolution vector field is 



where and N'^ denote the time derivative of N and A^" with respect to respectively, 
introduced in subsection 12.2. 1[ Note that while the normal-normal and normal-tangential 
parts of the spacetime Killing operator depend on M and M", its spatial projection does 
not. It is well defined even on a single spacelike hypersurface. Thus, by (|3.17l) . the 
integrand of the n-dimensional integral on the right-hand side of fl3.14l) is a°''^{N'V(aKb) — 



On the other hand, contrary to the bulk terms, the boundary integral (which we call 
the Poisson boundary term) apparently contains the canonical momentum explicitly, and 
not only through the various parts of the symmetric energy momentum tensor. However, 
if we take into account any of the boundary conditions (I3.12p and (I3.13p . then the Poisson 
boundary term can be rewritten purely in terms of the energy-momentum tensor. In fact, 
if i Km holds (both for {N,N^) and (iV,iV")), then by iVE^ = -N%Da^)U = -N^Ja 
the first and the third terms in the boundary integral of (13.141) cancel each other, and 
there remains only VaCr^-^iNNb — NNb). Similarly, using the explicit form of the spatial 
stress aab, we can write a'^^NNb = NN^-fi + N^N^Jb - iViVJ". Then, however, it is 
straightforward to check that the integrand of the boundary integral is an expression of 
the energy-momentum tensor, and the final form of (|3.14p is 



N + N^DaM - M^'DaN, 

N" + MD^N - ND^M - [M, N] 

NXah + D^aNb), 



a 



(3.15) 
(3.16) 
(3.17) 



NV(aKb)). 




= H ^N^'DaN - N^DaN, ND^N - ND^N - [N, N] "J + 
+ / (NW^aKb)-NW^aKb))a'''d^+ (f K'^^SacT'bK'dS. 



(3.18) 
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Using the boundary condition it is easy to check that the boundary integral is anti- 
symmetric in K'^K'^, as it should be because every other term in (13.181) changes sign if 
we interchange {N,N°-) and (iV, iV"). Similarly, if (13.131) holds, then the last term in the 
boundary integral in (l31iD is vanishing, and VaCr'^'^NNb = VaN'^Nfi + VaN'^N^Jb-NNv^Jb 
holds. Using this expression, (13.141) again takes the form (13.18p . 

Recalling that in the spacetime picture the Hamiltonian H[N, N^] is the flux integral 
of the Lorentz-covariant current T"-^Kb, it is natural to ask for the spacetime-covariant 
form of (13.18p . too. To derive this, let us decompose the Lie bracket [K,K]"- of the 
spacetime vector flelds K"^ and K"^ with respect to the foliation and evolution vector fleld 
above. It is 

[K,K]\ = tH'(^NV^bK,) - NV^bK,)^ + N-D,N - N^aN, (3.19) 

[K, K] V," = 2h-H' (iVV(feir,) - iVV(5irc)) - 

- {ND^N - ND^N) + [N, N] (3.20) 



Thus the new lapse and shift both in ^2.11\} and 13.141 ) appear as the lapse and shift parts 



of the Lie bracket of the spacetime vector fields up to the spacetime Killing operators. 
Substituting these into (13.18^ and using the notation H[K^] = H[N,N^], we obtain 

{h[K'] , H [K'] } = -h[ [K, W]+ {K,V(aKb) - KcV^aKb)) T^'dS + 

+ ^ jf {K'^K' - K-K') ^SacT^dS. (3.21) 

Therefore, the quasi-local Hamiltonians of the real scalar field do not form a Poisson 
algebra. The two roots of this failure are the non-Killing nature of the vector fields K"^ 
and K'' and the boundary integral. The latter is precisely the boundary integral of (13.21) . 
Our aim is to recover the whole of ( 13. 2p . 

To do this, let us calculate the time derivative of H[N, N"-] with respect to in the 
spacetime. It is 

^i7[iV,iV»]= J (^^R + J^R- + ^Na^'>hab)d^ + 

= J [liN+JaN'' + Na''\MXab + DiaMb)))d^ + 

+ ^H[M,M''],H[N,N^]^, (3.22) 

where in the first step we used the functional differentiability of if [TV, TV"], and in the 
second the canonical equations of motion with the Hamiltonian H[M, M°]. Finally, by the 
expression (13.18p of the Poisson bracket and the projections ( 13.15p -( l3TT7ll of the spacetime 
Killing operator, we obtain 

j^H [N, N^] = J CtcT'''V^aKb)dE + jf r ^BacT^'KbdS. (3.23) 
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Thus we recovered (13.21) . whose boundary integral appeared here as the Poisson boundary 
term. 



4 The Poisson boundary terms in GR 

4.1 The Poisson boundary terms 
4.1.1 The formal Poisson brackets 

Though in subsection 12.1.11 we defined the Poisson bracket for functionally dijferentiable 
functions on the phase space, we can define the formal Poisson bracket of any two con- 
straint functions C[N, N'^] and C[N, iV^] by the integral of their formal functional deriva- 
tives, independently of their functional differentiability. A lengthy but straightforward 
calculation gives (or see [HI ITT] ) 



C [N, N"] , C [N, N'']j = C ^N^DeN - N^'DeN , ND^N - ND^N - [N, N] ' 
- j 2vaP''' (^NDtN - ND,N - [N, N] dS - 



- (b 2p'''(VeN'DaNb - VeN'DaNb)dS 



i{l(i.-2A + ,^ 



XabX 



-v''{DaN){DtN'') + [D'N) {D,Nay}dS. 



(4.1) 



A well-known highly non-trivial property of the constraint functions is that in their formal 
Poisson bracket, the genuine n-dimensional integral is also a constraint function (with 
the new lapse NWeN - N^D^N and the new shift ND'^N - ND'^N - [N.N]"), and 
the remaining terms are all boundary integrals. It might also be interesting to note that 
the first two terms on the right together is just the basic Hamiltonian of subsection 12.3.21 
parameterized by the new lapse and shift. 



4.1.2 The main idea 

We learnt in subsection 13.2.21 that the Poisson boundary term in {H[^"-], H[K"-]} de- 
scribes, at least for appropriately chosen Killing fields the infinitesimal fiux of energy- 
momentum and angular momentum fiows through the hypersurface that is generated by 
Lie dragging S along the integral curves of However, if we could expect that the 
Poisson boundary term has a physical meaning in general relativity too, then (in ad- 
dition to the requirement of the functional differentiability of the Hamiltonian and the 
compatibility of the boundary conditions and the evolution equations) we would have a 
further condition that we could use to find the 'ultimate' Hamiltonian boundary term 
and the boundary conditions. 

Thus suppose for a moment that we already have the 'ultimate' Hamiltonian boundary 
term and the boundary conditions both for the canonical variables and the lapse and 
the shift; and hence H[N,N"-]^ given by p. 31) . is functionally differentiable. Then the 
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Poisson bracket of two such Hamiltonians is precisely the formal Poisson bracket of the 
two constraints with the same lapses and shifts, which is already given explicitly by 
f l4.ip . However, if our expectation is correct, then this Poisson bracket must be the sum 
of the 'correct' Hamiltonian (parameterized by the new lapse and shift) and another 
physical quantity (being analogous to the energy-momentum and angular momentum 
fluxes). Consequently, since both are gauge invariant, the Poisson boundary term must 
also be gauge invariant. In particular, this Poisson boundary term 

1. should depend only on N\s, N\s, N°-\s and N°'\s, but not on their normal derivatives, 

e.g. on v^DeN\s or v^DeN"-\s:, 

2. should depend only on the geometry of S, but not on the geometry of S at its 

boundary S; 

3. must be 2 + (n — l)-covariant and, in particular, it must be independent of the actual 

choice for the normals {t"", v"") to S. 

Therefore, we must check whether or not the boundary terms in (14. ip satisfy these criteria. 
This will be done by rewriting it in a form adapted to S, using the ideas and notions 
summarized in subsection I2.2.2[ 



4.1.3 The covariant form of the Poisson boundary term 



Using the definitions of subsection 12.2.21 by a systematic decomposition of every tensor 
field and derivative operator according to P^" = 11^ — v'^Vb we rewrite the integrand of 
every boundary integral in (14. ip . Our ultimate aim is to obtain a 2 + {n — l)-covariant 
form and, in particular, in terms of the spacetime vector fields K'^ := Nf^ + and 
K'^ := Nt^ + rather than the individual lapses and shifts. 

First, by a tedious but straightforward computation for the integrand Ii of the first 
boundary integral in (14. ip . we obtain 



2kJi 



+ 



+ 
+ 
+ 



K'5eK'' - K'5eK''j {rva - Aa) + 



VeN^'^DaNb - VeN%^DaNb) [A' - Tv'). 



(4.2) 



(To reproduce this formula (and the similar ones below), it seems useful to calculate 
first the projections of various quantities, e.g. 2KVaVbP"''^ 



XabV"-v'' + X 



2KVaP'^''qbc = Ac or 2Kp'"^qcaqdb = Tab - Qabir - XcdV^'v'^). Note that Xabv"-v^ cannot be 
expressed by quantities defined only on S.) 

Similarly, the integrand I2 of the second boundary integral is 
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+ (veK'{KK'') - v,K'{KK''))xcdv''v^ + 

+ (veK^'^DaNb - VeK'v^DaNb) {tv' - A') . (4.3) 

Before rewriting the integrand of the third boundary integral let us observe that the 
Ricci tensor and the curvature scalar appear in ly, just like in the 'constraint parts' v"-v^Gab 
and v'^Gab^l of the Einstein tensor of the n-dimensional intrinsic (spatial) geometry 
(S, hab)- Expressing these in terms of the curvature scalar and the extrinsic curvature 
Vab of S (and its derivative SeUab) , a direct but quite lengthy calculation gives that 

+ l{^R+ (i^afei^"' -ly'- TabT"' + t''))^ S .^K'^ + 

+ v''{DaN) (^A^K^ - tN^ - v''{DaN) (^A^K^ - riv) + 
+ (2A„A« - TXabV^v^Ys^K'^K'' + 
+ {K\5aKb) - K^{5aKb)yv + 

+ - iyq"^)t,(K'^5aK^ - K'^SaKb) + 

+ {5,Ka){S'Kb) + 

+ A''(^{6bK,yv,K'^ - {6bKa)v%,K'') + 

+ - {5bKa)t''t,K'^^ + 

+ A^Tcatb [K'^K^ - K^K^) . (4.4) 

Adding the three integrands and forming total ^a-divergences, the resulting expression 
can be written in the form 



21 



+ \{^R + ra^r'^' -r^- ^ab^"' + J^') ^e.^K^K'' + 

+ [AbK") {v^DeN + VeN'xcdv'v'' - A,N' - ^e'^K^ - 

+ 5a[{K''K'-K''K'')A,+ {p'^'-pq''')t''{K,K,-K,K,)^. (4.5) 

Since the third line is just the contraction of the curvature of the connection 5e in the 
normal bundle and the vector fields K"^ and K"^, moreover the fourth line is proportional 
to the trace of (12 .71) : one might attempt to rewrite (14.51) in a form containing the curvature 
F'^bcd of Ag. We show that this can indeed be done. 

To get terms such as the right hand side of ( 12.81) and (12.91) . let us rewrite the first as 
well as the fifth lines of ( 14. 5p as total (5e-divergences and terms with the derivative of the 
extrinsic curvature tensor. Re-expressing the second line in terms of the Ae-derivative 
operator we have 

2K{h + I2 + h) = -A ^ScdK^^K" - [A^K^) [AS') ^Sab + 

+ l{^R+ l^aby""' -V''- TabT^' + T^) ^S.dK'K'' + 

+ i (6aA, - 6,Aa + r.Veft - r^Ve,) [K^K' - K'^K') + 
+ [s^Q\e - 5,Q\e) ^e'diK'K'' - K'K'') + 
+ {A,K'>) (v'D.N + v.N'Xcdv'v'' - A.N'^'^ - 
- {A,K') (v'D.N + VeN^.dv'v'' - A.N'''^ + 

+ 5, - K^K') (A + ^e,,Qr) + [K^'K, - K'K,) ^e^/Q"'^ + 

+ _ j,q-^)t-{K,Kb - KcKb)). (4.6) 

Comparing its second, third and fourth lines with (I2.7p - (l2.10p . we find that these can be 
rewritten as 

= lK'K'-F'\f^E„S:lS, (4.7) 

by means of which we arrive at our final expression for the formal Poisson bracket of two 
constraint functions: 
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C[N,N^] ,C[N,N''] 
= C 



+ 



-{A^K^) (vW,N + v.N^Xcdv'v^ - N'A^i }d5. (4.8) 

Since the curvature F°'bcd is the pull back to S of the spacetime curvature 2-form ^R"-f,cd, 
in the physically important special case n = 3 the curvature term reduces to ^K"'K^eabcd 
R'^'^ef^^^ , and hence gives a Penrose-type charge integral of the spacetime curvature |29j . 
In a GHP spin frame (o^, l^) adapted to S this takes the form K'^ K'^e a' b' {"^ abcdo'^ — 
^ABC'D'O^ + ^{oa^b + i-aOb)) + C.C.. Thus, in particular, it is only the ^2 and 
the \l/3, but not the \l/o and \l/4 Weyl spinor components that are involved in the Poisson 
boundary term. 

4.1.4 Boundary conditions from the gauge invariance of the Poisson bound- 
ary term 

Clearly, the first three terms in the boundary integral of (14. 8p are manifestly 2 + (n — 1)- 
covariant; they depend only on the geometry of S and the value of the vector fields on 
S (but independent of the way in which they are extended off the boundary) , and they 
are invariant with respect to the change of the actual normals (t",f°) of iS. On the other 
hand, the last two lines contain 'bad' terms. Thus, we can ensure the gauge invariance of 
the Poisson boundary term (in the sense discussed in subsection 14.1.21) if we require the 
vanishing of the Ae-divergence of the vector fields and K"^. Obviously, AaK'^ = is a 
2+{n-l) covariant condition, and from A^AT" = SaK'^+Q'^abK^ = Nr-VaN^u+SaiU^N'') 
it is clear that it has infinitely many solutions on S: the condition AaK"- = specifies 
e.g. only the lapse in terms of the still completely freely specifiable shift N"-. 

To see the meaning of this condition, let us rewrite this into the form g'^^n^II^ 
V(c-Arf) = 0. This is one of the |n(n — 1) projected parts of the |(n+ l)(n + 2) spacetime 
Killing equations. Thus, AaK"- = is a weakening of the familiar spacetime Killing 
equations. Clearly, if the lapse part of K"^ is vanishing and the shift part N"- is tangent 
to S on S, then AaK"' = reduces to the condition SaN"" = discussed in subsection 
12X21 

To clarify its compatibility with the evolution equations, let us rewrite the canonical 
equation of motion for the metric hab in the spacetime. By (|3.17p . its right-hand side is 
just the projection to S of the spacetime Killing operator: 

hab = 2NXab + L^hab = 2P^P,^\/ (cKd) . 

Hence, the contraction of its restriction to S with the metric g"^ gives q°'^hab = 2AbK^. 
However, the left-hand side is proportional to the time derivative of the induced volume 
element on S\ £ei...e„„i = |g"''ga6£ei...e„_i = y'^''hab£ei...e„^i- Therefore, AaK"" = is 
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precisely the condition that the induced volume {n — l)-form on S is constant during the 
evolution. Thus the boundary condition (5eei...e„_i = for the configuration variables, 
found in the special cases and discussed in subsections 12.3.11 and 12.3.21 appears naturally 
in the general case, too. 



4.1.5 On Einstein-scalar systems 

The quasi-local phase space of the coupled Einstein-scalar system is the cotangent bundle 
T*Q(E) of the configuration space Q(S), the latter being the set of the pairs {hat, $), and 
endowed with the natural symplectic structure. The constraint of the coupled system is 
^[A^, A^"] := Ce[N, N"] + Hs[N, N""] = 0, where now ^ is denoted by Ce[N, N''] and 
Hs[N, N"-] is given by (13.51) . and the Hamiltonian is the sum of the Hamiltonians of the 
gravitational and the scalar sectors: H[N, N"-] = He[N, N'^] + Hs[N, iV"]. However, note 
that Hs[N,N"'] depends on the metric hat, which is now a configuration variable. Thus, 
assuming that both A'" := Nt°- + N°- and K"^ := Nt"- + N"- are Ag-divergence-free, the 
formal Poisson bracket of the Hamiltonians H[N, N"-] and H[N, iV"] is 



H[N,N''] , H[N,N''] 



Since F'^bcd 
as 



M Da 



■ C 



(4.9) 



N^D^N - N^D^N , ND^N - ND^N - [A^, N] ' 



R^bef^t^d^ the last three terms of the boundary integral can be written 



2 ' Vn — 1 n[n — 1) / 

where '^Cabcd, ^ Rah and ^ R are the spacetime Weyl and Ricci tensors and the cur- 
vature scalar, respectively, and the second line is proportional to the expression whose 
vanishing is just the Einstein equation. Thus, in particular for = 3 and 'on shell', the 
Poisson bracket of two Hamiltonians is the boundary integral of {/S.'^K'^){/S.eK^) '^^ah — 
iK'^K'^^^Cabcd^e'"^ + l^RK^K^-^Eab- Hence, the trace-free part of the spacetime Ricci 
tensor does not appear even in the presence of a scalar field. 

4.2 Quasi-local quantities from Poisson boundary terms? 

In subsection 14.1.21 we raised the idea that the Poisson boundary term should be the 
sum of the Hamiltonian boundary term (parameterized by the new lapse and shift) and 
terms analogous to the fiux of energy-momentum/angular momentum of matter fields. 



24 



and both must be gauge invariant. In the present subsection, we decompose the (gauge- 
invariant) Poisson boundary term in such a way that the boundary term of the 'improved' 
basic Hamiltonian (12.131) emerges naturally, even in a (slightly modified) gauge-invariant 
form. However, as we will see, in its gauge-invariant form it does not seem to yield a 
representation of the 'composition rule' of how the new lapse and shift are built from the 
old ones. 

We start with (14. 5p . and let us observe first that in its first line the derivative operator 
5e can be replaced by Ag; moreover, the first term in the second line can be written as 
\{5e6^K°-)-^eabK^ — \{^e^'^K"-)'^^abK^ up to a total (5e-divergence. Again, by forming 
total (5(j-divergences, the third line is written as {K^/S.f,K"- — K'^AeK^)Ah plus extrinsic 
curvature terms. Thus, for vector fields satisfying AaK"- = AaK'^ = 0, we obtain 



^H[N,N''] ,H[N,N'']^ = (4.10) 
= C In^D^N - N^DeN , ND^N - ND^N - [N, N] " 

+i(...^ir-2(^.X.)g--2ir^Q..Q-- 

Thus the first line in the boundary integral is just the boundary term of Hi[K°-] given 
by (I2T3D . in which K" is replaced by the 'commutator' K^'AeK" - K^'A^K"-. Since, 
assuming AaK°- = AaK"' = 0, the integral of (|4.51) is invariant with respect to the change 
of the basis (t°',v"') in the normal bundle of S ('S'0(1,1) boost-gauge invariance'), the 
decomposition of the integrand of (14.10p to the Hamiltonian boundary term (the first line) 
and to the rest can be made in a boost-gauge-invariant way, too. In fact, by the Hodge 
decomposition (see e.g. [30]) the connection 1-form Af> is the sum of an exact, a co-exact 
and a harmonic 1-form on iS: A^ = Sga + ag + uig, respectively, and this decomposition 
is unique. Here the function a : 5 — > R is unique up to an additive constant, ae is 
(5e-divergence free, while Ug is both 5e-divergence free and satisfies 6[aUJb] = 0. (The first 
represents the pure gauge part of Ae, the co-exact part yields curvature, and the harmonic 
part only holonomy, but no curvature.) Thus by 6[aAh] = 6[aCib] it is only the co-exact 
part of Ae that appears in (14.51) . and hence we can substitute A^ by Oe in (14.101) . too, 
yielding a manifestly boost gauge invariant form of the 'improved' basic Hamiltonian. 
Another interpretation of the above Hodge decomposition is that it provides a 'natural' 
gauge fixing, using only the intrinsic geometry of S. 

Unfortunately, however, this Hamiltonian boundary term does not seem to represent 
the 'composition rule' of the lapses and shifts in a correct way. Indeed, the lapse-shift 
parts of K'^AeK" - K^AgK'' are not the new lapse N'^DgN - N'^DgN and the new shift 
ND^N - ND^N - [N, N]" that appear in the constraint function. 
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5 Conclusions 



We learnt from the quasi-local canonical formulation of the scalar field that the Poisson 
boundary terms represent energy-momentum and angular momentum fluxes out from 
and into the localized system, i.e. they are well defined physical quantities. Here we raise 
the idea that the same may be expected in general relativity, too, and hence the Poisson 
boundary term in GR must be gauge invariant in every sense. We showed that this 
requirement yields the condition for the lapse and shift that the spacetime vector field 
that they determine must be divergence free with respect to a Sen-type connection on the 
boundary. This condition is a part of the spacetime Killing equations. This yields that 
the evolution equations preserve the volume form induced on the boundary. Therefore, 
keeping the induced volume form fixed as a condition seems to be the part of the (yet 
unknown) 'ultimate' boundary conditions for the canonical variables. This implies that 
the quasi-local constraint algebra that we found earlier (and discussed in subsection l2.3.ip 
is probably the 'correct' one, completing the point (iii) of subsection I2.1.T1 

We also found arguments both in favour of and against the Hamiltonian boundary 
term in Hi[K"-]. It appears naturally as a part of the Poisson boundary term and its 
50(1,1) boost-gauge dependence can be cured, but, without additional restrictions on 
the lapse and the shift it does not yield functionally differentiable Hamiltonian, and it 
does not represent the composition law for the lapses and shifts in a correct way. 

The author is grateful to Robert Beig, Jorg Frauendiener, Helmut Friedrich and Niall 
O Murchadha for the useful remarks and discussions on the canonical formulation of gen- 
eral relativity and the role of boundary conditions, to Edward Anderson and the referees 
for their bibliographic remarks, as well as the Erwin Schrodinger Institute, Vienna, for 
hospitality. This work was partially supported by the Hungarian Scientific Research Fund 
(OTKA) grant K67790. 
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